Method for adaptively adjusting control parameters of BBMC-based speed regulation system of asynchronous motor

ABSTRACT

A mathematical model between the optimization objects and the optimization objectives is established, with finite-time control parameters as optimization objects, and a total harmonic distortion of an output voltage, a deviation signal of a capacitor voltage and a deviation signal of an output current of BBMC as optimization objectives. A multi-objective optimization satisfaction function and a multi-objective optimization fitness function are established. Multiple sets of optimal finite-time control parameters are iteratively optimized using an adaptive wolf pack algorithm. Functional relationships between respective optimal finite-time control parameters and the corresponding actual output current of BBMC are obtained using a numerical fitting method. Optimal control parameters corresponding to any load are obtained according to the obtained functional relationships. Therefore, the control parameters are adjusted in real time according to the functional relationships and the load of the BBMC-based speed regulation system or the actual output current of BBMC.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application is a continuation of International Application No. PCT/CN2019/098333 with a filling date of Jul. 30, 2019, which claims the benefit of priority from Chinese Patent Application No. 201910171144.6 with a filing date of Mar. 7, 2019. The content of the aforementioned applications, including any intervening amendments thereto, are incorporated herein by reference.

TECHNICAL FIELD

The present application relates to power electronics, and more particularly to a method for adaptively adjusting control parameters of a BBMC-based speed regulation system of an asynchronous motor.

BACKGROUND

As a novel power converter, a Buck-Boost matrix converter (BBMC) has a high voltage transfer ratio and is capable of directly outputting high-quality sine waves, which is applicable for speed regulation systems of asynchronous motors. However, research shows that, in the BBMC-based speed regulation system of the asynchronous motor, control parameters of the speed regulation system of the asynchronous motor are required to be optimized as the load of the BBMC-based speed regulation system changes, so that indicators such as the total harmonic distortion of the output voltage in the BBMC achieve an optimal collaboration, thereby achieving a high-performance speed regulation for the speed regulation system.

Therefore, it is of great significance to determine optimal control parameters according to the magnitude of the load (i.e, magnitude of an actual output current of the BBMC), and further study a change law between the optimal control parameters and the actual output current of the BBMC, so that the relevant control parameters are adjusted in real time according to the load changes, thereby realizing the high-performance speed regulation.

SUMMARY

In order to solve the above-mentioned technical problems, the present disclosure provides a method for adaptively adjusting control parameters of a BBMC-based speed regulation system of an asynchronous motor with a high precision and good operation effect.

The technical solutions of the disclosure are described as follows.

The present disclosure provides a method for adaptively adjusting control parameters of a BBMC-based speed regulation system of an asynchronous motor, comprising:

(1) taking finite-time control parameters as optimization objects, a total harmonic distortion of an output voltage of a BBMC, a deviation signal of a capacitor voltage of the BBMC and a deviation signal of an output current of the BBMC as optimization objectives, establishing a mathematical model between the optimization objects and the optimization objectives for the BBMC-based speed regulation system of the asynchronous motor;

(2) selecting a load to allow an actual output current of the BBMC to be a certain value; and establishing a multi-objective optimization satisfaction function and a multi-objective optimization fitness function for the BBMC-based speed regulation system;

(3) iteratively optimizing the finite-time control parameters using an adaptive wolf pack optimization algorithm, so that an optimal collaboration among the total harmonic distortion of the output voltage of the BBMC, the deviation signal of the capacitor voltage and the deviation signal of the output current is obtained, thereby obtaining a set of optimal finite-time control parameters; adjusting the actual output current of the BBMC and repeating the step (2) and the step (3) to obtain n sets of optimal finite-time control parameters; and

(4) obtaining functional relationships between the optimal finite-time control parameters and the corresponding actual output current of the BBMC using a numerical fitting method according to the obtained n sets of optimal finite-time control parameters and the corresponding actual output currents of the BBMC; and determining the optimal finite-time control parameters corresponding to any load of the BBMC-based speed regulation system according to the obtained functional relationships.

In an embodiment, the mathematical model between the optimization objects and the optimization objectives is established through steps of:

(1.1) taking a capacitor voltage u_(C), an inductor current i_(L) and an output current i_(l) as control variables of the BBMC-based speed regulation system, establishing a state differential function of the BBMC-based speed regulation system as follows:

$\begin{matrix} \left\{ {\begin{matrix} {{\overset{.}{i}}_{L} = {{- \frac{u_{C}}{L}} + {\frac{u_{C} + u_{D}}{L}d}}} \\ {{\overset{.}{u}}_{C} = {\frac{i_{L}}{C} - \frac{i_{1}}{C} - {\frac{i_{L}}{C}d}}} \\ {{\overset{.}{i}}_{1} = {\frac{i_{1}R}{L_{1}} - \frac{u_{C}}{L_{1}} + \frac{u_{DZ}}{L_{1}}}} \end{matrix};} \right. & (1) \end{matrix}$

wherein u_(C) is the capacitor voltage; i_(L) is the inductor current; i_(l) is the output current of the BBMC; u_(D) is the voltage in a DC side of the BBMC; u_(DZ) is the voltage in a common terminal of the single-phase stator windings of the asynchronous motor; L and C are respectively the inductance parameter and the capacitance parameter at an inverter side of the BBMC; R and L₁ are respectively the equivalent resistance and the equivalent inductance of the single-phase stator windings of the asynchronous motor; and d is the duty cycle of a power switch tube in the BBMC;

(1.2) obtaining an equation of the duty cycle d of the power switch tube in the BBMC according to the equation (1) and a finite-time control algorithm:

$\begin{matrix} {{d = \frac{\begin{matrix} {{{CL}_{1}{u_{0}\left( {{L\left( {{- {K_{1}\left( {{{sat}_{\alpha_{1}}\left( \lambda_{1} \right)} + \lambda_{1}} \right)}} - {K_{2}\left( {{{sat}_{\alpha_{2}}\left( \lambda_{2} \right)} + \lambda_{2}} \right)}} \right)} + {u_{C}u_{D}}} \right)}} +} \\ {{L_{1}{Li}_{1}^{2}u_{C}} + {{{CLu}_{D}\left( {u_{C} + u_{D}} \right)}\left( {{Ri}_{1} - u_{C} + u_{DZ}} \right)}} \end{matrix}}{{{CL}_{1}\left( {u_{C} + u_{D}} \right)}\left( {u_{D}^{2} + i_{L}^{2}} \right)}};} & (2) \end{matrix}$

wherein sat is the saturation function; λ₁=(Li_(L) ²+C(u_(C)+u_(D))²−Li_(L) _(ref) ²−C(u_(C) _(ref) +u_(D))²)/2; λ₂=u_(D)i_(L)−(u_(C)=u_(D))i_(l)−u_(D)i_(Lref)+(u_(Cref)+u_(D))i_(lref); u_(Cref) is the reference value of the capacitor voltage; i_(Lref) is the reference value of the inductor current; i_(lref) is the reference value of the output current in the BBMC; K₁, K₂, α₁ and α₂ are the finite-time control parameters; and α₂=2α₁/(1+α₁);

(1.3) respectively obtaining analytical equations of the output voltage u and the output current i_(l) according to the equation (1):

$\begin{matrix} {{u = {\left( {{\exp\left( {{\ln\sqrt{\frac{{R^{2}C} - {4L_{1}}}{4L_{1}{LC}}}} + {\frac{R}{2}\sqrt{\frac{{R^{2}C} + {4L_{1}}}{4L_{1}^{2}C}}}} \right)} + {\left( {1 - d} \right){Ri}_{L}}} \right)\left( {{\cos\left( \sqrt{\frac{{R^{2}C} + {4L_{1}}}{4{RL}_{1}{LC}}t} \right)} + {\sqrt{\frac{{R^{2}C} + {4L_{i}}}{L_{1}^{2}C}}{\sin\left( {\sqrt{\frac{{R^{2}C} + {4L_{1}}}{4{RL}_{1}{LC}}}t} \right)}}} \right)}},} & (3) \\ {i_{1} = {\left( {\exp\left( {{\ln\left( {2{C^{2}\left( {1 - d} \right)}i_{L}} \right)} + {\frac{1}{2C}\sqrt{\frac{4}{L_{1}C} + \frac{R^{2}}{L_{1}^{2}}}} + \frac{R}{L_{1}}} \right)} \right)\left( {{{\frac{R}{L_{1}}{\cos\left( {\sqrt{\frac{4}{L_{1}C} + \frac{R^{2}}{L_{1}^{2}}}t} \right)}} + {\sin\left( {\sqrt{\frac{4}{L_{1}C} + \frac{R^{2}}{L_{1}^{2}}}t} \right)}};} \right.}} & (4) \end{matrix}$

(1.4) obtaining the total harmonic distortion THD of the output voltage u according to a definition of the total harmonic distortion:

$\begin{matrix} {{{THD} = \frac{\sqrt{\begin{matrix} {\left( {{\exp\left( {{TH} + G} \right)} - {2\;\omega\;{KRi}_{L}d} + {\left( {R + {L_{1}H}} \right)/\left( {8L_{1}\omega^{2}} \right)}} \right)^{2} +} \\ \left( {{\exp\left( {{TH} + G} \right)} - {3\;\omega\;{KRi}_{L}d} - {\left( {R + {L_{1}H}} \right)/\left( {16L_{1}\omega^{2}} \right)}} \right)^{2} \end{matrix}}}{{\exp\left( {{TH} + G} \right)} - {\omega\;{KRi}_{L}d} - {\left( {R + {L_{1}H}} \right)/\left( {2L_{1}\omega^{2}} \right)}}};} & (5) \end{matrix}$

wherein

${H = \sqrt{\frac{{R^{2}C} + {4L_{1}}}{L_{1}^{2}C}}},\mspace{14mu}{K = \sqrt{\frac{{R^{2}C} - {4L_{1}}}{4L_{1}{LC}}}},$ G=ln(R/2), T is the period of the output voltage of the BBMC; and ω is the angular frequency of the output voltage; and

(1.5) respectively obtaining a deviation Δu_(c) between the capacitor voltage u_(c) and an ideal voltage value U_(e) under a load and a deviation Δi between the output current i_(l) and an ideal current I_(e) under a load according to the analytical equations of the output voltage u and the output current i_(l) obtain ed from the equation (3) and the equation (4):

$\begin{matrix} {{{\Delta\; u_{c}} = {{u_{c} - U_{e}} = {{\left( {{\exp\left( {{\ln\sqrt{\frac{{R^{2}C} - {4L_{1}}}{4L_{1}{LC}}}} + \frac{R}{2} - \sqrt{\frac{{R^{2}C} + {4L_{1}}}{4L_{1}^{2}C}}} \right)} + {\left( {1 - d} \right){Ri}_{L}}} \right)\left( {{\cos\left( {\sqrt{\frac{{R^{2}C} + {4L_{1}}}{4{RL}_{1}{LC}}}t} \right)} + {\sqrt{\frac{{R^{2}C} + {4L_{i}}}{L_{1}^{2}C}}{\sin\left( {\sqrt{\frac{{R^{2}C} + {4L_{1}}}{4{RL}_{1}{LC}}}t} \right)}}} \right)} + U_{DZ} - {Ue}}}},} & (6) \\ {{\Delta\; i} = {i_{1} = {I_{e} = {\left( {\exp\left( {{\ln\left( {2{C^{2}\left( {1 - d} \right)}i_{L}} \right)} + {\frac{1}{2C}\sqrt{\frac{4}{L_{1}C} + \frac{R^{2}}{L_{1}^{2}}}} + \frac{R}{L_{1}}} \right)} \right)\left( {{\frac{R}{L_{1}}{\cos\left( {\sqrt{\frac{4}{L_{1}C} + \frac{R^{2}}{L_{1}^{2}}}t} \right)}} + {\sin\left( {\sqrt{\frac{4}{L_{1}C} + \frac{R^{2}}{L_{1}^{2}}}t} \right)} - {I_{e}.}} \right.}}}} & (7) \end{matrix}$

In an embodiment, the multi-objective optimization satisfaction function is established through steps of:

(2.1.1) respectively establishing satisfaction functions of optimization objectives THD, Δu_(C) and Δi, wherein

a satisfaction function ƒ₁ of an optimization objective THD is expressed as shown in equation (8):

$\begin{matrix} {f_{1}\left\{ {\begin{matrix} 1 & {{THD} \leq {THD}^{\prime}} \\ \frac{c_{1}}{\sqrt{\left( {{THD} - {THD}^{\prime}} \right)}} & {{THD} > {THD}^{\prime}} \end{matrix};} \right.} & (8) \end{matrix}$

a satisfaction function ƒ₂ of an optimization objective Δu_(C) is expressed as shown in equation (9):

$\begin{matrix} {f_{2} = \left\{ {\begin{matrix} 1 & {{\Delta\; u_{C}} \leq {\Delta\; u_{C}^{\prime}}} \\ \frac{c_{2}}{\sqrt{\left( {{\Delta\; u_{C}} - {\Delta\; u_{C}^{\prime}}} \right)}} & {{\Delta\; u_{C}} > {\Delta\; u_{C}^{\prime}}} \end{matrix};} \right.} & (9) \end{matrix}$

a satisfaction functions ƒ₃ of an optimization objective Δi is expressed as shown in equation (10):

$\begin{matrix} {f_{3} = \left\{ {\begin{matrix} 1 & {{\Delta\; i} \leq {\Delta\; i^{\prime}}} \\ \frac{c_{3}}{\sqrt{\left( {{\Delta\; i} - {\Delta\; i^{\prime}}} \right)}} & {{\Delta\; i} > {\Delta\; i^{\prime}}} \end{matrix};} \right.} & (10) \end{matrix}$

wherein THD′, Δu_(C)′ and Δi′ are respectively thresholds of the optimization objectives THD, Δu_(C) and Δi; c₁, c₂ and c₃ are coefficients of satisfaction, wherein c₁>0, c₂>0 and c₃>0; and

(2.1.2) establishing the multi-objective optimization satisfaction function ƒ of the optimization objectives THD, Δu_(C) and Δi as shown in equation (11): ƒ=k ₁ƒ₁ +k ₂ƒ₂ +k ₃ƒ₃  (11);

wherein k, k₂ and k₃ are respectively weight coefficients of the optimization objectives THD, Δu_(C) and Δi; and k₁+k₂+k₃=1; and the multi-objective optimization fitness function is established through steps of:

(2.2.1) comparing satisfaction values of the optimization objectives with corresponding satisfaction thresholds; if a satisfaction value ƒ_(j) (j=1, 2, 3) of each of the optimization objectives is smaller than a corresponding satisfaction threshold M_(j) (j=1, 2, 3), the satisfaction value ƒ_(j) (j=1, 2, 3) is multiplied by a dynamic penalty factor b_(j), wherein the satisfaction thresholds are respectively

${M_{1} = {\frac{{THD} - {THD}^{\prime}}{THD}}},\mspace{14mu}{M_{2\;} = {{\frac{{\Delta\; u_{C\;}} - {\Delta\; u_{C}^{\prime}}}{\Delta\; u_{C}}}\mspace{14mu}{and}}}$ ${M_{3} = {\frac{{\Delta\; i} - {\Delta\; i^{\prime}}}{\Delta\; i}}};$ and the corresponding dynamic penalty factors are respectively

${b_{1} = {\log_{\frac{1}{2}}\left( {\frac{{THD} - {THD}^{\prime}}{THD}} \right)}},\mspace{14mu}{b_{2} = {{\log_{\frac{1}{2}}\left( {\frac{{\Delta\; u_{C\;}} - {\Delta\; u_{C}^{\prime}}}{\Delta\; u_{C}}}\mspace{11mu} \right)}\mspace{20mu}{and}}}$ ${b_{3} = {\log_{\frac{1}{2}}\left( {\frac{{\Delta\; i} - {\Delta\; i^{\prime}}}{\Delta\; i}} \right)}};$

otherwise, if the satisfaction value ƒ_(j) (j=1, 2, 3) of each of the optimization objectives is greater than or equal to the corresponding satisfaction threshold M_(j) (j=1, 2, 3), setting the dynamic penalty factor b_(j) as 1;

(2.2.2) establishing the multi-objective optimization fitness function ƒ, as shown in equation (12): ƒ_(s) =k ₁ b ₁ƒ₁ +k ₂ b ₂ƒ₂ +k ₃ b ₃ƒ₃  (12).

In an embodiment, the finite-time control parameters in the BBMC are iteratively optimized using the adaptive wolf pack optimization algorithm through steps of:

(3.1) taking the selected actual output current of the BBMC as a judgment reference value of the adaptive wolf pack optimization algorithm;

(3.2) initializing parameters, wherein the parameters comprise: population size N representing N sets of the finite-time control parameters, the maximum number k_(max) of iterations and the multi-objective optimization fitness function ƒ_(s) representing a concentration S(i) of prey's smell of a wolf pack;

(3.3) configuring a random direction and a random distance of wolves in the wolf pack and obtaining a position Xi (K₁, K₂, α₁) of the i-th wolf;

(3.4) obtaining the concentration of prey's smell of the wolf pack to obtain a set of the optimization objects; and obtaining the corresponding THD, Δu_(C) and Δi according to the optimization objects;

(3.5) identifying a wolf with the highest concentration of prey's smell in the wolf pack as an optimal wolf, and retaining the concentration of prey's smell and position X_(m)(K₁, K₂, α₁) of the optimal wolf;

(3.6) eliminating N/10 wolves in the wolf pack with a small concentration of prey's smell, and randomly generating the same number of new wolves in a solution space to update the wolf pack;

(3.7) determining whether the maximum number of iterations is reached; if yes, outputting the position X_(m) (K₁, K₂, α₁) of the optimal wolf to output an optimal solution of the finite-time control parameters K₁, K₂ and α₁; otherwise, adding 1 to the number of iterations, and returning to the step (3.3);

(3.8) determining whether n sets of optimal finite-time control parameters are obtained; if no, changing the actual output current of the BBMC at a certain interval, and returning to the step (3.1); and

(3.9) outputting the n sets of optimal finite-time control parameters and the corresponding actual output currents of the BBMC.

In an embodiment, in the step (4), the numerical fitting method is a least square method; the functional relationships comprise a functional relationship between an optimal finite-time control parameter K₁ and the actual output current i of the BBMC, a functional relationship between an optimal finite-time control parameter K₂ and the actual output current i of the BBMC, and a functional relationship between an optimal finite-time control parameter α₁ and the actual output current i of the BBMC;

wherein the functional relationship between the optimal finite-time control parameter K₁ and the actual output current i of the BBMC is expressed as shown in equation (13): ƒ_(K1)(i)=A ₁ i ⁴ +A ₂ i ³ +A ₃ i ² +A ₄ i+A ₅  (13);

the functional relationship between the optimal finite-time control parameter K₂ and the actual output current i of the BBMC is expressed as shown in equation (14): ƒ_(K2)(i)=B ₁ i ⁴ +B ₂ i ³ +B ₃ i ² +B ₄ i+B ₅  (14);

the functional relationship between the optimal finite-time control parameter α₁ and the actual output current i of the BBMC is expressed as shown in equation (15): ƒ_(α1)(i)=C ₁ exp(−((i−C ₂)/C ₃)²)+C ₄ exp(−((i−C ₅)/C ₆)²)+C ₇ exp(−((i−C ₈)/C ₉)²)+C ₁₀ exp(−((i−C ₁₁)²)+C ₁₃ exp(−((i−C ₁₄)/C ₁₄)/C ₁₅)²)  (15);

wherein ƒ_(K1)(i), ƒ_(K2)(i) and ƒ_(α1)(i) are respectively functions of the optimal finite-time control parameters K₁, K₂ and α₁; A₁, A₂, A₃, A₄ and A₅ are respectively coefficients of ƒ_(K1)(i); B₁, B₂, B₃, B₄ and B₅ are respectively coefficients of ƒ_(K2)(i); C₁, C₂, C₃, C₄, C₅, C₆, C₇, C₈, C₉, C₁₀, C₁₁, C₁₂, C₁₃, C₁₄ and C₁₅ are respectively coefficients of ƒ_(α1)(i); and the coefficients are obtained using a Matlab™ analysis software according to the least square method; and

the optimal finite-time control parameters corresponding to any load of the BBMC-based speed regulation system are determined according to the obtained functional relationships.

The disclosure has the following beneficial effects.

The present disclosure provides a method for adaptively adjusting control parameters of a BBMC-based speed regulation system of an asynchronous motor, in which the BBMC-based speed regulation system of the asynchronous motor is adopted as a frequency converter, which is controlled by a finite-time control algorithm. Firstly, taking finite-time control parameters as optimization objects, and a total harmonic distortion THD of an output voltage, the deviation signal Δu_(c) of the capacitor voltage and the deviation signal Δi of output current as optimization objectives, a mathematical model between the optimization objects and the optimization objectives is established. Then, a load is selected to allow an actual output current of the BBMC to be a certain value, and a multi-objective optimization satisfaction function and a multi-objective optimization fitness function are established; a set of optimal finite-time control parameters are obtained using an adaptive wolf pack algorithm, and the actual output current of the BBMC is adjusted at a certain interval to obtain n sets of optimal finite-time control parameters. Finally, functional relationships between the optimal finite-tithe control parameters and the corresponding actual output current of the BBMC are obtained using a numerical fitting method. The control parameters are adjusted in real time according to the actual load of the speed regulation system, so that the speed regulation system achieves an optimal operation effect under any load.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 shows a topology of a BBMC-based speed regulation system according to an embodiment of the present disclosure.

FIG. 2 is a flowchart of a method for adaptively adjusting control parameters of the BBMC-based speed regulation system of an asynchronous motor according to the present disclosure.

FIG. 3 is a flowchart of an adaptive wolf pack optimization algorithm according to the present disclosure.

FIGS. 4A-C show curves of fitting functions of optimal finite-time control parameters according to the present disclosure.

DETAILED DESCRIPTION OF EMBODIMENTS

The present disclosure will be further described in detail below with reference to the accompanying drawings and embodiments.

As shown in FIG. 1, provided is a topology of a speed regulation system according to the present disclosure, where a BBMC is adopted as a frequency converter, and a three-phase asynchronous motor is adopted as a load motor; the BBMC includes a rectifier side and an inverter side, where the rectifier side is a three-phase PWM rectifier circuit for rectifying three-phase AC voltage into PWM-modulated DC voltage; and the inverter side is a three-phase Buck-Boost inverter, which consists of three Buck-Boost DC/DC converters with the same structure.

As shown in FIG. 2, provided is a flowchart of a method for adaptively adjusting control parameters of the BBMC-based speed regulation system of the asynchronous motor, including the following steps:

(1) taking finite-time control parameters as optimization objects, and a total harmonic distortion of an output voltage of the BBMC, a deviation signal of a capacitor voltage of the BBMC and a deviation signal of an output current of the BBMC as optimization objectives, a mathematical model between the optimization objects and the optimization objectives for the BBMC-based speed regulation system of the asynchronous motor is established;

where the mathematical model between the optimization objects and the optimization objectives is established through steps of:

(1.1) taking a capacitor voltage u_(c), an inductor current i_(L) and an output current i_(l) as control variables of the speed regulation system, a state differential function of the speed regulation system is established as follows:

$\begin{matrix} \left\{ {\begin{matrix} {{\overset{.}{i}}_{L} = {{- \frac{u_{C}}{L}} + {\frac{u_{C} + u_{D}}{L}d}}} \\ {{\overset{.}{u}}_{C} = {\frac{i_{L}}{C} - \frac{i_{1}}{C} - {\frac{i_{L}}{C}d}}} \\ {{\overset{.}{i}}_{1} = {\frac{i_{1}R}{L_{1}} - \frac{u_{C}}{L_{1}} + \frac{u_{DZ}}{1}}} \end{matrix};} \right. & (1) \end{matrix}$

where u_(C) is the capacitor voltage; i_(L) is the inductor current; i_(l) is the output current of the BBMC; u_(D) is the voltage in a DC side of the BBMC; u_(DZ) is the voltage in a common terminal of single-phase stator windings of the asynchronous motor; L and C are respectively the inductance parameter and the capacitance parameter at an inverter side of the BBMC; R and L₁ are respectively the equivalent resistance and the equivalent inductance of the single-phase stator windings of the asynchronous motor; and d is the duty cycle of a power switch tube in the BBMC;

(1.2) an equation of the duty cycle d of the power switch tube in the BBMC is obtained according to the equation (1) and a finite-time control algorithm:

$\begin{matrix} {{d = \frac{\begin{matrix} {{{CL}_{1}{u_{0}\left( {{L\left( {{- {K_{1}\left( {{{sat}_{\alpha_{1}}\left( \lambda_{1} \right)} + \lambda_{1}} \right)}} - {K_{2}\left( {{{sat}_{\alpha_{2}}\left( \lambda_{2} \right)} + \lambda_{2}} \right)}} \right)} + {u_{C}u_{D}}} \right)}} +} \\ {{L_{1}{Li}_{1}^{2}u_{C}} + {{{CLu}_{D}\left( {u_{C} + u_{D}} \right)}\left( {{Ri}_{1} - u_{C} + u_{DZ}} \right)}} \end{matrix}}{{{CL}_{1}\left( {u_{C} + u_{D}} \right)}\left( {u_{D}^{2} + i_{L}^{2}} \right)}};} & (2) \end{matrix}$

where sat is the saturation function; λ₁=(Li_(L) ²+C(u_(C)+u_(D))²−Li_(L) _(ref) ²−C(u_(C) _(ref) +u_(D))²)/2; λ₂=u_(D)i_(L)−(u_(C)+u_(D))i_(l)−u_(D)i_(Lref)+(u_(Cref)+u_(D))i_(lref); u_(Cref) is the reference value of the capacitor voltage; i_(Lref) is the reference value of the inductor current; i_(lref) is the reference value of the output current in the BBMC; K₁, K₂, α₁ and α₂ are the finite-time control parameters; and α₂=2α₁/(1+α₁);

(1.3) analytical equations of the output voltage u and the output current i_(l) are respectively obtained according to the equation (1):

$\begin{matrix} {{u = {\left( {{\exp\left( {{\ln\sqrt{\frac{{R^{2}C} - {4L_{1}}}{4L_{1}{LC}}}} + {\frac{R}{2}\sqrt{\frac{{R^{2}C} + {4L_{1}}}{4L_{1}^{2}C}}}} \right)} + {\left( {1 - d} \right){Ri}_{L}}} \right)\left( {{\cos\left( {\sqrt{\frac{{R^{2}C} + {4L_{1}}}{4{RL}_{1}{LC}}}t} \right)} + {\sqrt{\frac{{R^{2}C} + {4L_{i}}}{L_{1}^{2}C}}\sin\mspace{11mu}\left( {\sqrt{\frac{{R^{2}C} + {4L_{1}}}{4{RL}_{1}{LC}}}t} \right)}} \right)}},} & (3) \\ {{i_{1} = {\left( {\exp\left( {{\ln\left( {2{C^{2}\left( {1 - d} \right)}i_{L}} \right)} + {\frac{1}{2C}\sqrt{\frac{4}{L_{1}C} + \frac{R^{2}}{L_{1}^{2}}}} + \frac{R}{L_{1}}} \right)} \right)\left( {{\frac{R}{L_{1}}{\cos\left( {\sqrt{\frac{4}{L_{1}C} + \frac{R^{2}}{L_{1}^{2}}}t} \right)}} + {\sin\left( {\sqrt{\frac{4}{L_{1}C} + \frac{R^{2}}{L_{1}^{2}}}t} \right)}} \right)}};} & (4) \end{matrix}$

(1.4) the total harmonic distortion THD of the output voltage u is obtained according to a definition of the total harmonic distortion:

$\begin{matrix} {{{THD} = \frac{\sqrt{\begin{matrix} {\left( {{\exp\left( {{TH} + G} \right)} - {2\;\omega\;{KRi}_{L}d} + {\left( {R + {L_{1}H}} \right)/\left( {8L_{1}\omega^{2}} \right)}} \right)^{2} +} \\ \left( {{\exp\left( {{TH} + G} \right)} - {3\;\omega\;{KRi}_{L}d} - {\left( {R + {L_{1}H}} \right)/\left( {16L_{1}\omega^{2}} \right)}} \right)^{2} \end{matrix}}}{{\exp\left( {{TH} + G} \right)} - {\omega\;{KRi}_{L}d} - {\left( {R + {L_{1}H}} \right)/\left( {2L_{1}\omega^{2}} \right)}}},} & (5) \end{matrix}$

where

${H = \sqrt{\frac{{R^{2}C} + {4L_{1}}}{L_{1}^{2}C}}},\mspace{14mu}{K = \sqrt{\frac{{R^{2}C} - {4L_{1}}}{4L_{1}{LC}}}},$ G=ln(R/2), T is the period of the output voltage of the BBMC; and ω is the angular frequency of the output voltage; and

(1.5) a deviation Δu_(c) between the capacitor voltage u_(c) and an ideal voltage value U_(e) under a certain load and a deviation Δi between the output current i_(l) and an ideal current I_(e) under a certain load are respectively obtained according to the analytical equations of the output voltage u and the output current i_(l) obtained from the equation (3) and the equation (4):

$\begin{matrix} {{{\Delta\; u_{c}} = {{u_{c} - U_{e}} = {{\left( {{\exp\left( {{\ln\sqrt{\frac{{R^{2}C} - {4L_{1}}}{4L_{1}{LC}}}} + \frac{R}{2} - \sqrt{\frac{{R^{2}C} + {4L_{1}}}{4L_{1}^{2}C}}} \right)} + {\left( {1 - d} \right){Ri}_{L}}} \right)\left( {{\cos\left( {\sqrt{\frac{{R^{2}C} + {4L_{1}}}{4{RL}_{1}{LC}}}t} \right)} + {\sqrt{\frac{{R^{2}C} + {4L_{i}}}{L_{1}^{2}C}}{\sin\left( {\sqrt{\frac{{R^{2}C} + {4L_{1}}}{4{RL}_{1}{LC}}}t} \right)}}} \right)} + U_{DZ} - {Ue}}}},} & (6) \\ {{{\Delta\; i} = {i_{1} = {I_{e} = {{\left( {\exp\left( {{\ln\left( {2{C^{2}\left( {1 - d} \right)}i_{L}} \right)} + {\frac{1}{2C}\sqrt{\frac{4}{L_{1}C} + \frac{R^{2}}{L_{1}^{2}}}} + \frac{R}{L_{1}}} \right)} \right)\left( {{\frac{R}{L_{1}}{\cos\left( {\sqrt{\frac{4}{L_{1}C} + \frac{R^{2}}{L_{1}^{2}}}t} \right)}} + {\sin\left( {\sqrt{\frac{4}{L_{1}C} + \frac{R^{2}}{L_{1}^{2}}}t} \right)}} \right)} - I_{e}}}}};} & (7) \end{matrix}$

(2) a load with a certain magnitude is selected to allow an actual output current of the BBMC to be a certain value; and a multi-objective optimization satisfaction function and a multi-objective optimization fitness function for the speed regulation system are established;

(2.1) the multi-objective optimization satisfaction function is established through the following steps:

(2.1.1) satisfaction functions of optimization objectives THD, Δu_(C) and Δi are respectively established, where

a satisfaction function ƒ₁ of the optimization objective THD is expressed as shown in equation (8):

$\begin{matrix} {f_{1}\left\{ {\begin{matrix} 1 & {{THD} \leq {THD}^{\prime}} \\ \frac{c_{1}}{\sqrt{\left( {{THD} - {THD}^{\prime}} \right)}} & {{THD} > {THD}^{\prime}} \end{matrix};} \right.} & (8) \end{matrix}$

a satisfaction function ƒ₂ of the optimization objective Δu_(C) is expressed as shown in equation (9):

$\begin{matrix} {f_{2} = \left\{ {\begin{matrix} 1 & {{\Delta u_{C}} \leq {\Delta\; u_{C}^{\prime}}} \\ \frac{c_{2}}{\sqrt{\left( {{\Delta\; u_{C}} - {\Delta\; u_{C}^{\prime}}} \right)}} & {{\Delta u_{C}} > {\Delta\; u_{C}^{\prime}}} \end{matrix};} \right.} & (9) \end{matrix}$

a satisfaction functions ƒ₃ of the optimization objective Δi is expressed as shown in equation (10):

$\begin{matrix} {f_{3} = \left\{ {\begin{matrix} 1 & {{\Delta\; i} \leq {\Delta\; i^{\prime}}} \\ \frac{c_{2}}{\sqrt{\left( {{\Delta\; i} - {\Delta\; i^{\prime}}} \right)}} & {{\Delta\; i} > {\Delta\; i^{\prime}}} \end{matrix};} \right.} & (10) \end{matrix}$

where THD′, Δu_(C)′ and Δi′ are respectively thresholds of the optimization objectives THD, Δu_(C) and Δi; c₁, c₂ and c₃ are coefficients of satisfaction, wherein c₁>0, c₂>0 and c₃>0; and

(2.1.2) the multi-objective optimization satisfaction function ƒ of the optimization objectives THD, Δu_(C) and Δi is established as shown in equation (11): ƒ=k ₁ƒ₁ +k ₂ƒ₂ +k ₃ƒ₃  (11);

where k, k₂ and k₃ are respectively weight coefficients of the optimization objectives THD, Δu_(C) and Δi; and k₁+k₂+k₃=1;

(2.2) the multi-objective optimization fitness function is established through steps of:

(2.2.1) a satisfaction value of each of the optimization objectives is compared with its corresponding satisfaction threshold; if a satisfaction value ƒ_(j) (j=1, 2, 3) of each of the optimization objectives is smaller than its corresponding satisfaction threshold M_(j) (j=1, 2, 3), the satisfaction value ƒ_(j) (j=1, 2, 3) of each of the optimization objectives is multiplied by a corresponding dynamic penalty factor b_(j), where the satisfaction thresholds are respectively

${M_{1} = {\frac{{THD} - {THD}^{\prime}}{THD}}},{{M_{2} = {{{\frac{{\Delta u_{C}} - {\Delta\; u_{C}^{\prime}}}{\Delta u_{C}}}\mspace{14mu}{and}\mspace{14mu} M_{3}} = {\frac{{\Delta\; i} - {\Delta\; i^{\prime}}}{\Delta\; i}}}};}$ and the corresponding dynamic penalty factors are respectively

${{b_{1} = {\log_{\frac{1}{2}}\left( {\frac{{THD} - {THD}^{\prime}}{THD}} \right)}},{b_{2} = {{\log_{\frac{1}{2}}\left( {\frac{{\Delta u_{C}} - {\Delta\; u_{C}^{\prime}}}{\Delta u_{C}}} \right)}\mspace{14mu}{and}}}}\;$ ${b_{3} = {\log_{\frac{1}{2}}\left( {\frac{{\Delta\; i} - {\Delta\; i^{\prime}}}{\Delta\; i}} \right)}};$

otherwise, i.e., if the satisfaction value ƒ_(j) (j=1, 2, 3) of each of the optimization objectives is greater than or equal to its corresponding satisfaction threshold M_(j) (j=1, 2, 3), the corresponding dynamic penalty factor b_(j) is set as 1;

(2.2.2) the multi-objective optimization fitness function ƒ₅ is established as shown in equation (12): ƒ_(s) =k ₁ b ₁ƒ₁ +k ₂ b ₂ƒ₂ +k ₃ b ₃ƒ₃  (12);

(3) the finite-time control parameters are iteratively optimized using an adaptive wolf pack optimization algorithm, so that an optimal collaboration among the total harmonic distortion of the output voltage, the deviation signal of the capacitor voltage and the deviation signal of the output current is obtained, thereby obtaining a set of optimal finite-time control parameters; the actual output current of the BBMC is adjusted and the step (2) and the step (3) are repeated to obtain n sets of optimal finite-time control parameters.

Referring to FIG. 3, the finite-time control parameters in the BBMC are iteratively optimized by adopting the adaptive wolf pack optimization algorithm through the following steps:

(3.1) the selected actual output current of the BBMC is taken as a judgment reference value of the adaptive wolf pack optimization algorithm;

(3.2) parameters are initialized, where the parameters include: population size N representing N sets of the finite-time control parameters, the maximum number k_(max) of iterations and the multi-objective optimization fitness function ƒ_(s) representing a concentration S(i) of prey's smell of a wolf pack;

(3.3) a random direction and a random distance of wolves in the wolf pack are configured and a position Xi (K₁, K₂, α₁) of the i-th wolf is obtain;

(3.4) the concentration of prey's smell of the wolf pack is obtained to obtain a set of the optimization objects; and the corresponding THD, Δu_(C) and Δi are obtained according to the optimization objects;

(3.5) a wolf with the highest concentration of prey's smell in the wolf pack is identified as an optimal wolf, and the concentration of prey's smell and position X_(m) (K₁, K₂, α₁) of the optimal wolf are retained;

(3.6) N/10 wolves in the wolf pack with a small concentration of prey's smell are eliminated, and the same number of new wolves are randomly generated in a solution space to update the wolf pack;

(3.7) whether the maximum number of iterations is reached is determined; if yes, the position X_(m) (K₁, K₂, α₁) of the optimal wolf is outputted, i.e., an optimal solution of the finite-time control parameters K₁, K₂ and α₁ are outputted; otherwise, 1 is added to the number of iterations, and the process is returned to step (3.3);

(3.8) whether n sets of optimal finite-time control parameters are obtained is determined; if no, the actual output current of the BBMC is changed at a certain interval, and the process is returned to the step (3.1); and

(3.9) the n sets of optimal finite-time control parameters and the corresponding actual output currents of the BBMC are outputted.

(4) Functional relationships between the optimal finite time control parameters and the actual output current of the BBMC are obtained using a numerical fitting method according to the obtained n sets of optimal finite-time control parameters and the corresponding actual output currents of the BBMC; the numerical fitting method is a least square method; the functional relationships include a functional relationship between an optimal finite-time control parameter K₁ and the actual output current i of the BBMC, a functional relationship between an optimal finite-time control parameter K₂ and the actual output current i of the BBMC, and a functional relationship between an optimal finite-time control parameter α₁ and the actual output current i of the BBMC;

where the functional relationship between the optimal finite-time control parameters K₁ and the actual output current i of the BBMC is expressed as shown in equation (13): ƒ_(K1)(i)=A ₁ i ⁴ +A ₂ i ³ +A ₃ i ² +A ₄ i+A ₅  (13);

the functional relationship between the optimal finite-time control parameters K₂ and the actual output current i of the BBMC is expressed as shown in equation (14): ƒ_(K2)(i)=B ₁ i ⁴ +B ₂ i ³ +B ₃ i ² +B ₄ i+B ₅  (14);

the functional relationship between the optimal finite-time control parameters α₁ and the actual output current i of the BBMC is expressed as shown in equation (15): ƒ_(α1)(i)=C ₁ exp(−((i−C ₂)/C ₃)²)+C ₄ exp(−((i−C ₅)/C ₆)²)+C ₇ exp(−((i−C ₈)/C ₉)²)+C ₁₀ exp(−((i−C ₁₁)²)+C ₁₃ exp(−((i−C ₁₄)/C ₁₄)/C ₁₅)²)  (15);

where ƒ_(K1)(i), ƒ_(K2)(i) and ƒ_(α1)(i) are respectively functions of the optimal finite-time control parameters K₁, K₂ and α₁; A₁, A₂, A₃, A₄ and A₅ are respectively coefficients of ƒ_(K1)(i); B₁, B₂, B₃, B₄ and B₅ are respectively coefficients of ƒ_(K2)(i); C₁, C₂, C₃, C₄, C₅, C₆, C₇, C₈, C₉, C₁₀, C₁₁, C₁₂, C₁₃, C₁₄ and C₁₅ are respectively coefficients of ƒ_(α1)(i); and the coefficients are obtained using a Matlab™ analysis software according to the least square method; and

the optimal finite-time control parameters corresponding to any load of the speed regulation system are determined according to the obtained functional relationships.

FIGS. 4A-C show curves of fitting functions of optimal finite-time control parameters according to the present disclosure. A rated power, a rated voltage, and a rated current of the speed regulation system are set as follows: P_(N)=15 kW, U_(N)=380 V, I_(N)=30.1 A. 30 sets of load data are taken at a certain interval within a range of the rated load of the speed regulation system, i.e., 30 output current data are selected at a certain interval within the range of the rated output current of the BBMC. For example, 30 output current data are selected by taking 1.4 A as an initial value and increasing with an interval of 0.4 A, and the optimal finite-time control parameters corresponding to each output current are calculated using the adaptive wolf pack algorithm, as as shown in Table 1:

TABLE 1 Output currents of the BBMC and the corresponding optimal finite-time control parameters Current/A 1.4 1.8 2.2 2.6 3.0 K₁ 0.05331 0.07200 0.09193 0.11297 0.13500 K₂ 0.03554 0.04800 0.11288 0.25313 0.40003 α₁ 0.12500 0.14285 0.16667 0.14285 0.12500 Current/A 3.4 3.8 4.2 4.6 5.0 K₁ 0.15794 0.18169 0.20619 0.23136 0.25715 K₂ 0.55294 0.71130 0.87460 1.04241 1.21435 α₁ 0.11111 0.08333 0.07692 0.07142 0.07692 Current/A 5.4 5.8 6.2 6.6 7.0 K₁ 0.28351 0.31039 0.33775 0.36555 0.39376 K₂ 1.39007 1.56927 1.75167 1.93702 2.12511 α₁ 0.08333 0.09090 0.16667 0.16667 0.14285 Current/A 7.4 7.8 8.2 8.6 9.0 K₁ 0.42235 0.44216 0.46770 0.48816 0.50852 K₂ 2.31573 2.49593 2.61804 2.75446 2.89015 α₁ 0.11111 0.09090 0.08653 0.08142 0.07692 Current/A 9.4 9.8 10.2 10.6 11.0 K₁ 0.52877 0.54892 0.56899 0.58897 0.60887 K₂ 3.02515 3.15952 3.29329 3.42650 3.55919 α₁ 0.07893 0.08111 0.11500 0.12285 0.14256 Current/A 11.4 11.8 12.2 12.6 13.0 K₁ 0.62870 0.64846 0.66815 0.68778 0.70735 K₂ 3.69137 3.82309 3.95436 4.08522 4.21568 α₁ 0.16667 0.12500 0.10000 0.08142 0.07692

According to the 30 sets of output currents of the BBMC and the corresponding optimal finite-time control parameters obtained from Table 1, the corresponding functional relationships are obtained using the numerical fitting method; where the numerical fitting method is the least square method; the functional relationships include a functional relationship between the optimal finite-time control parameter K₁ and the actual output current i of the BBMC, a functional relationship between the optimal finite-time control parameter K₂ and the actual output current i of the BBMC, and a functional relationship between the optimal finite-time control parameter α₁ and the actual output current i of the BBMC.

(a) The functional relationship between the optimal finite-time control parameter K₁ and the actual output current i of the BBMC is expressed as shown in equation (13): ƒ_(K1)(i)=A ₁ i ⁴ +A ₂ i ³ +A ₃ i ² +A ₄ i+A ₅  (13);

where A₁=4.163×10⁻⁵, A₂=−0.001382, A₃=0.01484, A₄=0.00125 and A₅=0.02925.

(b) The functional relationship between the optimal finite-time control parameter K₂ and the actual output current i of the BBMC is expressed as shown in equation (14): ƒ_(K2)(i)=B ₁ i ⁴ +B ₂ i ³ +B ₃ i ² +B ₄ i+B ₅  (14);

where B₁=0.0004296, B₂=−0.01408, B₃=0.1537, B₄=−0.2423 and B₅=0.0802.

(c) The functional relationship between the optimal finite-time control parameter α₁ and the actual output current i of the BBMC is expressed as shown in equation (15): ƒ_(α1)(i)=C ₁ exp(−((i−C ₂)/C ₃)²)+C ₄ exp(−((i−C ₅)/C ₆)²)+C ₇ exp(−((i−C ₈)/C ₉)²)+C ₁₀ exp(−((i−C ₁₁)²)+C ₁₃ exp(−((i−C ₁₄)/C ₁₄)/C ₁₅)²)  (15);

where C₁=0.148, C₂=11.29, C₃=1.84, C₄=0.07538, C₅6.413, C₆=0.4914, C₇=0.07578, C₈=2.162, C₉=1.123, C₁₀=0.08371, C₁₁=2.678, C₁₂=4.058, C₁₃=0.09212, C₁₄=7.344 and C₁₅=1.663.

The coefficients of the functional relationships corresponding to the optimal finite-time control parameters are obtained using a Matlab™ analysis software according to the least square method; based on the functional relationship between the optimal control parameters and the actual output current of the BBMC, the control parameters are adjusted in real time according to the actual load of the speed regulation system), so that the speed regulation system can achieve an optimal operation effect under any load.

Described above is only a preferred embodiment of the present disclosure. It should be noted that any improvement and variation made by those skilled in the art without departing from the spirit of the present disclosure shall fall within the scope of the appended claims. 

What is claimed is:
 1. A method for adaptively adjusting control parameters of a Buck-Boost matrix converter (BBMC) based speed regulation system of an asynchronous motor, comprising: (1) taking finite-time control parameters as optimization objects, a total harmonic distortion of an output voltage of a BBMC, a deviation signal of a capacitor voltage of the BBMC and a deviation signal of an output current of the BBMC as optimization objectives, establishing a mathematical model between the optimization objects and the optimization objectives for the speed BBMC-based regulation system of the asynchronous motor; wherein the mathematical model between the optimization objects and the optimization objectives is established through steps of: (1.1) taking a capacitor voltage u_(C), an inductor current i_(L) and an output current i_(l) as control variables of the BBMC-based speed regulation system, establishing a state differential function of the BBMC-based speed regulation system as follows: $\begin{matrix} \left\{ {\begin{matrix} {\overset{.}{i_{L}} = {\frac{u_{C}}{L} + {\frac{u_{C} + u_{D}}{L}d}}} \\ {\overset{.}{u_{C}} = {\frac{i_{L}}{C} - \frac{i_{1}}{C} - {\frac{i_{L}}{C}d}}} \\ {\overset{.}{i_{1}} = {\frac{i_{1}R}{L_{1}} - \frac{u_{C}}{L_{1}} + \frac{u_{DZ}}{L_{1}}}} \end{matrix};} \right. & (1) \end{matrix}$ wherein u_(C) is the capacitor voltage; i_(L) is the inductor current; i_(l) is the output current of the BBMC; u_(D) is the voltage in a DC side of the BBMC; u_(DZ) is the voltage in a common terminal of single-phase stator windings of the asynchronous motor; L and C are respectively the inductance parameter and the capacitance parameter at an inverter side of the BBMC; R and L₁ are respectively the equivalent resistance and the equivalent inductance of the single-phase stator windings of the asynchronous motor; and d is the duty cycle of a power switch tube in the BBMC; (1.2) obtaining an equation of the duty cycle d of the power switch tube in the BBMC according to the equation (1) and a finite-time control algorithm: $\begin{matrix} {{d = \frac{\left( {{{CL}_{1}{u_{D}\left( {{L\left( {{- {K_{1}\left( {{{sat}_{\alpha_{1}}\left( \lambda_{1} \right)} + \lambda_{1}} \right)}} - {K_{2}\left( {{{sat}_{\alpha_{2}}\left( \lambda_{2} \right)} + \lambda_{2}} \right)}} \right)} + {u_{C}u_{D}}} \right)}} + {L_{1}L_{1}^{2}u_{C}} + {{{CLu}_{D}\left( {u_{C} + u_{D}} \right)}\left( {{Ri}_{1} - u_{C} + u_{DZ}} \right)}} \right.}{\left( {{{CL}_{1}\left( {u_{C} + u_{D}} \right)}\left( {u_{D}^{2} + i_{L}^{2}} \right)} \right)}};} & (2) \end{matrix}$ wherein sat is the saturation function; λ₁=(Li_(L) ²+C(u_(C)+u_(D))²−Li_(L) _(ref) ²−C(u_(C) _(ref) +u_(D))²)/2; λ₂=u_(D)i_(L)−(u_(C)=u_(D))i_(l)−u_(D)i_(Lref)+(u_(Cref)+u_(D))i_(lref); u_(Cref) is the reference value of the capacitor voltage; i_(Lref) is the reference value of the inductor current; i_(lref) is the reference value of the output current in the BBMC; K₁, K₂, α₁ and α₂ are the finite-time control parameters; and α₂=2α₁/(1+α₁); (1.3) respectively obtaining analytical equations of the output voltage u and the output current i_(l) according to the equation (1): $\begin{matrix} {{u = {\left( {{\exp\left( {{\ln\sqrt{\frac{{R^{2}C} - {4\; L_{1}}}{4L_{1}{LC}}}} + {\frac{R}{2}\sqrt{\frac{{R^{2}C} + {4\; L_{1}}}{L_{1}^{2}C}}}} \right)} + {\left( {1 - d} \right){Ri}_{L}}} \right)\left( {{\cos\left( {\sqrt{\frac{{R^{2}C} + {4\; L_{1}}}{4{RL}_{1}{LC}}}t} \right)} + {\sqrt{\frac{{R^{2}C} + {4\; L_{1}}}{L_{1}^{2}C}}{\sin\left( {\sqrt{\frac{{R^{2}C} + {4\; L_{1}}}{4{RL}_{1}{LC}}}t} \right)}}} \right)}},} & (3) \\ {{i_{1} = {\left( {\exp\left( {{\ln\left( {2{C^{2}\left( {1 - d} \right)}i_{L}} \right)} + {\frac{1}{2\; C}\sqrt{\frac{4}{L_{1}C} + \frac{R^{2}}{L_{1}^{2}}}} + \frac{R}{L_{1}}} \right)} \right)\left( {{\frac{R}{L_{1}}{\cos\left( {\sqrt{\frac{4}{L_{1}C} + \frac{R^{2}}{L_{1}^{2}}}t} \right)}} + {\sin\left( {\sqrt{\frac{4}{L_{1}C} + \frac{R^{2}}{L_{1}^{2}}}t} \right)}} \right)}};} & (4) \end{matrix}$ (1.4) obtaining the total harmonic distortion THD of the output voltage u according to a definition of the total harmonic distortion: $\begin{matrix} {{THD} = \frac{\sqrt{\begin{matrix} {\left( {{\exp\left( {{TH} + G} \right)} - {2\omega\;{KRi}_{L}d} + {\left( {R + {L_{1}H}} \right)\text{/}\left( {8L_{1}\omega^{2}} \right)}} \right)^{2} +} \\ \left( {{\exp\left( {{TH} + G} \right)} - {3\omega\;{KRi}_{L}d} - {\left( {R + {L_{1}H}} \right)\text{/}\left( {16\; L_{1}\omega^{2}} \right)}} \right)^{2} \end{matrix}}}{{\exp\left( {{TH} + G} \right)} - {\omega\;{KRi}_{L}d} - {\left( {R + {L_{1}H}} \right)\text{/}\left( {2L_{1}\omega^{2}} \right)}}} & (5) \end{matrix}$ wherein ${H = \sqrt{\frac{{R^{2}C} + {4\; L_{1}}}{L_{1}^{2}C}}},{K = \sqrt{\frac{{R^{2}C} - {4\; L_{1}}}{4L_{1}{LC}}}},$ G=ln(R/2), T is the period of the output voltage of the BBMC; and ω is the angular frequency of the output voltage; and (1.5) respectively obtaining a deviation Δu_(c) between the capacitor voltage u_(c) and an ideal voltage value U_(e) under a load and a deviation Δi between the output current i_(l) and an ideal current I_(e) under a load according to the analytical equations of the output voltage u and the output current i_(l) obtained from the equation (3) and the equation (4): $\begin{matrix} {{{\Delta u_{c}} = {u_{c} - {{U_{e}\left( {{\exp\left( {{\ln\sqrt{\frac{{R^{2}C} - {4\; L_{1}}}{4L_{1}{LC}}}} + \frac{R}{2} - \sqrt{\frac{{R^{2}C} + {4\; L_{1}}}{L_{1}^{2}C}}} \right)} + {\left( {1 - d} \right){Ri}_{L}}} \right)}\left( {{\cos\left( {\sqrt{\frac{{R^{2}C} + {4\; L_{1}}}{4{RL}_{1}{LC}}}t} \right)} + {\sqrt{\frac{{R^{2}C} + {4\; L_{1}}}{L_{1}^{2}C}}{\sin\left( {\sqrt{\frac{{R^{2}C} + {4\; L_{1}}}{4{RL}_{1}{LC}}}t} \right)}}} \right)} + U_{DZ} - U_{e}}},} & (6) \\ {{{\Delta\; i} = {{i_{1} - I_{e}} = {{\left( {\exp\left( {{\ln\left( {2{C^{2}\left( {1 - d} \right)}i_{L}} \right)} + {\frac{1}{2\; C}\sqrt{\frac{4}{L_{1}C} + \frac{R^{2}}{L_{1}^{2}}}} + \frac{R}{L_{1}}} \right)} \right)\left( {{\frac{R}{L_{1}}{\cos\left( {\sqrt{\frac{4}{L_{1}C} + \frac{R^{2}}{L_{1}^{2}}}t} \right)}} + {\sin\left( {\sqrt{\frac{4}{L_{1}C} + \frac{R^{2}}{L_{1}^{2}}}t} \right)}} \right)} - I_{e}}}};} & (7) \end{matrix}$ (2) arbitrarily selecting a load of the BBMC-based speed regulation system to allow an actual output current of the BBMC to be a value; and establishing, corresponding to the selected load and the actual output current of the BBMC, a multi-objective optimization satisfaction function and a multi-objective optimization fitness function for the BBMC-based speed regulation system; wherein the multi-objective optimization satisfaction function is established through steps of: (2.1.1) respectively establishing satisfaction functions of optimization objectives THD, Δu_(C) and Δi, wherein a satisfaction function ƒ₁ of an optimization objective THD is expressed as shown in equation (8): $\begin{matrix} {f_{1} = \left\{ {\begin{matrix} 1 & {{THD} \leq {THD}^{\prime}} \\ \frac{c_{1}}{\sqrt{\left( {{THD} - {THD}^{\prime}} \right)}} & {{THD} > {THD}^{\prime}} \end{matrix};} \right.} & (8) \end{matrix}$ a satisfaction function ƒ₂ of an optimization objective Δu_(C) is expressed as shown in equation (9): $\begin{matrix} {f_{2} = \left\{ {\begin{matrix} 1 & {{\Delta u_{C}} \leq {\Delta\; u_{C}^{\prime}}} \\ \frac{c_{2}}{\sqrt{\left( {{\Delta\; u_{C}} - {\Delta\; u_{C}^{\prime}}} \right)}} & {{\Delta u_{C}} > {\Delta\; u_{C}^{\prime}}} \end{matrix};} \right.} & (9) \end{matrix}$ a satisfaction functions ƒ₃ of an optimization objective Δi is expressed as shown in equation (10): $\begin{matrix} {f_{3} = \left\{ {\begin{matrix} 1 & {{\Delta\; i} \leq {\Delta\; i^{\prime}}} \\ \frac{c_{2}}{\sqrt{\left( {{\Delta\; i} - {\Delta\; i^{\prime}}} \right)}} & {{\Delta\; i} > {\Delta\; i^{\prime}}} \end{matrix};} \right.} & (10) \end{matrix}$ wherein THD′, Δu_(C)′ and Δi′ are respectively thresholds of the optimization objectives THD, Δu_(C) and Δi; c₁, c₂ and c₃ are coefficients of satisfaction, wherein c₁>0, c₂>0 and c₃>0; and (2.1.2) establishing the multi-objective optimization satisfaction function ƒ of the optimization objectives THD, Δu_(C) and Δi as shown in equation (11): ƒ=k ₁ƒ₁ +k ₂ƒ₂ +k ₃ƒ₃  (11); wherein k, k₂ and k₃ are respectively weight coefficients of the optimization objectives THD, Δu_(C) and Δi; and k₁+k₂+k₃=1; and the multi-objective optimization fitness function is established through steps of: (2.2.1) comparing satisfaction values of the optimization objectives with corresponding satisfaction thresholds; if a satisfaction value ƒ_(j) (j=1, 2, 3) of each of the optimization objectives is smaller than a corresponding satisfaction threshold M_(j) (j=1, 2, 3), the satisfaction value ƒ_(j) (j=1, 2, 3) is multiplied by a dynamic penalty factor b_(j), wherein the satisfaction thresholds are respectively ${M_{1} = {\frac{{THD} - {THD}^{\prime}}{THD}}},{{M_{2} = {{{\frac{{\Delta u_{C}} - {\Delta\; u_{C}^{\prime}}}{\Delta u_{C}}}\mspace{14mu}{and}\mspace{14mu} M_{3}} = {\frac{{\Delta\; i} - {\Delta\; i^{\prime}}}{\Delta\; i}}}};}$ and the corresponding dynamic penalty factors are respectively ${{b_{1} = {\log_{\frac{1}{2}}\left( {\frac{{THD} - {THD}^{\prime}}{THD}} \right)}},{b_{2} = {{\log_{\frac{1}{2}}\left( {\frac{{\Delta u_{C}} - {\Delta\; u_{C}^{\prime}}}{\Delta u_{C}}} \right)}\mspace{14mu}{and}}}}\;$ ${b_{3} = {\log_{\frac{1}{2}}\left( {\frac{{\Delta\; i} - {\Delta\; i^{\prime}}}{\Delta\; i}} \right)}};$ otherwise, if the satisfaction value ƒ_(j) (j=1, 2, 3) of each of the optimization objectives is greater than or equal to its corresponding satisfaction threshold M_(j) (j=1, 2, 3), setting the dynamic penalty factor b_(j) as 1; (2.2.2) establishing the multi-objective optimization fitness function ƒ_(s) as shown in equation (12): ƒ_(s) =k ₁ b ₁ƒ₁ +k ₂ b ₂ƒ₂ +k ₃ b ₃ƒ₃  (12); (3) iteratively optimizing, using an adaptive wolf pack optimization algorithm, finite-time control parameters calculated from the multi-objective optimization satisfaction function and the multi-objective optimization fitness function, so that an optimal collaboration among the total harmonic distortion of the output voltage of the BBMC, the deviation signal of the capacitor voltage and the deviation signal of the output current is obtained, thereby obtaining a set of optimal finite-time control parameters corresponding to the selected load and the actual output current of the BBMC; repeating the step (2) and the step (3) to obtain n sets of optimal finite-time control parameters corresponding to n loads and n actual output currents of the BBMC; wherein the finite-time control parameters in the BBMC are iteratively optimized using the adaptive wolf pack optimization algorithm through steps of: (3.1) taking the selected actual output current of the BBMC as a judgment reference value of the adaptive wolf pack optimization algorithm; (3.2) initializing parameters, wherein the parameters comprise: population size N representing N sets of the finite-time control parameters, the maximum number k_(max) of iterations and the multi-objective optimization fitness function ƒ_(s) representing a concentration S(i) of prey's smell of a wolf pack; (3.3) configuring a random direction and a random distance of wolves in the wolf pack, and obtaining a position Xi (K₁, K₂, α₁) of the i-th wolf; (3.4) obtaining the concentration of prey's smell of the wolf pack to obtain a set of the optimization objects; and obtaining the corresponding THD, Δu_(C) and Δi according to the optimization objects; (3.5) identifying a wolf with the highest concentration of prey's smell in the wolf pack as an optimal wolf, and retaining the concentration of prey's smell and position X_(m)(K₁, K₂, α₁) of the optimal wolf; (3.6) eliminating N/10 wolves in the wolf pack with a small concentration of prey's smell, and randomly generating the same number of new wolves in a solution space to update the wolf pack; (3.7) determining whether the maximum number of iterations is reached; if yes, outputting the position X_(m) (K₁, K₂, α₁) of the optimal wolf to output an optimal solution of the finite-time control parameters K₁, K₂ and α₁; otherwise, adding 1 to the number of iterations, and returning to the step (3.3); (3.8) determining whether n sets of optimal finite-time control parameters are obtained; if no, changing the actual output current of the BBMC at a certain interval, and returning to the step (3.1); and (3.9) outputting the n sets of optimal finite-time control parameters and the corresponding actual output currents of the BBMC; (4) obtaining functional relationships between the n sets of optimal finite-time control parameters and the corresponding n loads and n actual output currents of the BBMC using a numerical fitting method; (5) determining an actual load of the speed regulation system in operation and determining a corresponding actual output current of the BBMC based on the actual load; and (6) adjusting control parameters of the speed regulation system to be optimal finite-time control parameters which are calculated based on the functional relationships from the actual load of the speed regulation system and the corresponding actual output current of the BBMC.
 2. The method of claim 1, wherein in the step (4), the numerical fitting method is a least square method; the functional relationships comprise a functional relationship between an optimal finite-time control parameter K₁ and the actual output current i of the BBMC, a functional relationship between an optimal finite-time control parameter K₂ and the actual output current i of the BBMC, and a functional relationship between an optimal finite-time control parameter α₁ and the actual output current i of the BBMC; wherein the functional relationship between the optimal finite-time control parameter K₁ and the actual output current i of the BBMC is expressed as shown in equation (13): ƒ_(K1)(i)=A ₁ i ⁴ +A ₂ i ³ +A ₃ i ² +A ₄ i+A ₅  (13); the functional relationship between the optimal finite-time control parameter K₂ and the actual output current i of the BBMC is expressed as shown in equation (14): ƒ_(K2)(i)=B ₁ i ⁴ +B ₂ i ³ +B ₃ i ² +B ₄ i+B ₅  (14); the functional relationship between the optimal finite-time control parameter α₁ and the actual output current i of the BBMC is expressed as shown in equation (15): ƒ_(α1)(i)=C ₁ exp(−((i−C ₂)/C ₃)²)+C ₄ exp(−((i−C ₅)/C ₆)²)+C ₇ exp(−((i−C ₈)/C ₉)²)+C ₁₀ exp(−((i−C ₁₁)²)+C ₁₃ exp(−((i−C ₁₄)/C ₁₄)/C ₁₅)²)  (15); wherein ƒ_(K1)(i), ƒ_(K2)(i) and ƒ_(α1)(i) are respectively functions of the optimal finite-time control parameters K₁, K₂ and α₁; A₁, A₂, A₃, A₄ and A₅ are respectively coefficients of ƒ_(K1)(i); B₁, B₂, B₃, B₄ and B₅ are respectively coefficients of ƒ_(K2)(i); C₁, C₂, C₃, C₄, C₅, C₆, C₇, C₈, C₉, C₁₀, C₁₁, C₁₂, C₁₃, C₁₄ and C₁₅ are respectively coefficients of ƒ_(α1)(i); and the coefficients are obtained using Matlab™ analysis software according to the least square method; and the optimal finite-time control parameters corresponding to any load of the BBMC-based speed regulation system are determined according to the obtained functional relationships. 